
BIREFRINGENCE MEASUREMENT AT DEEP-ULTRAVIOLET 

WAVELENGTHS 



Technical Field 

5 This application relates to precise measurement of birefringence properties of 

optical elements, including optical elements that are components of systems that use 
deep ultraviolet (DUV) wavelengths. 

Background 

Many important optical materials exhibit birefringence. Birefringence means 

10 that different linear polarizations of light travel at different speeds through the 

material. These different polarizations are most often considered as two components 
of the polarized light, one being orthogonal to the other. 

Birefringence is an intrinsic property of many optical materials, and may also 
be induced by external forces. Retardation or retardance represents the integrated 

15 effect of birefringence acting along the path of a light beam traversing the sample. If 
the incident light beam is linearly polarized, two orthogonal components of the 
polarized light will exit the sample with a phase difference, called the retardance. The 
fundamental unit of retardance is length, such as nanometers (nm). It is frequently 
convenient, however, to express retardance in units of phase angle (waves, radians, or 

20 degrees), which is proportional to the retardance (nm) divided by the wavelength of 
the light (nm). An "average" birefringence for a sample is sometimes computed by 
dividing the measured retardation magnitude by the thickness of the sample. 

Oftentimes, the term "birefringence" is interchangeably used with and carries 
the same meaning as the term "retardance." Thus, unless stated otherwise, those 

25 terms are also interchangeably used below. 

The two orthogonal polarization components described above are parallel to 
two orthogonal axes, which are determined by the sample and are respectively called 
the "fast axis" and the "slow axis." The fast axis is the axis of the material that aligns 
with the faster moving component of the polarized light through the sample. 

30 Therefore, a complete description of the retardance of a sample along a given optical 



path requires specifying both the magnitude of the retardance and its relative angular 

orientation of the fast (or slow) axis of the sample. 

The need for precise measurement of birefringence properties has become 

increasingly important in a number of technical applications. For instance, it is 
5 important to specify linear birefringence (hence, the attendant induced retardance) in 

optical elements that are used in high-precision instruments employed in 

semiconductor and other industries. 

Moreover, the optical lithography industry is currently transitioning to the use 

of very short exposure wavelengths for the purpose of further reducing line weights 
10 (conductors, etc.) in integrated circuits, thereby to enhance performance of those 

circuits. In this regard, the next generation of optical lithography tools will use laser 

light having a wavelength of about 157 nanometers, which wavelength is often 

referred to as deep ultraviolet or DUV. 

It is important to precisely determine the retardance properties of optical 
15 elements or components that are used in systems, such as lithography tools, that 

employ DUV. Such a component may be, for example, a calcium fluoride (CaF 2 ) lens 

of a scanner or stepper. Since the retardance of such a component is a characteristic 

of both the component material as well as the wavelength of light penetrating the 

material, a system for measuring retardance properties must operate with a DUV light 
20 source and associated components for detecting and processing the associated light 

signals. 

The magnitude of the measured retardance of an optical element is a function 
of the thickness of the element, the thickness being measured in the direction that the 
light propagates through the sample. For example, a CaF2 optical element will have 
25 an intrinsic birefringence of about 12 nm for every centimeter (cm) of thickness. 
Consequently, for example, a 10 cm-thick CaF 2 element will have a relatively high 
birefringence level of about 120 nanometers, which is about three-quarters of a 157 
nm DUV wavelength. 



Summary of the Invention 

The present invention is directed to systems and methods for precisely 
measuring birefringence properties of optical elements, especially those elements that 
are used in DUV applications. 

In one preferred embodiment, the system includes two photoelastic modulators 
(PEM) located on opposite sides of the sample. Each PEM is operable for modulating 
the polarity of a light beam that passes though the sample. The system also includes a 
polarizer associated with one PEM, an analyzer associated with the other PEM, and a 
detector for measuring the intensity of the light after it passes through the PEMs, the 
polarizer, and the analyzer. 

As one aspect of the invention, embodiments and methods are described for 
determining birefringence properties across a wide range. For example, an 
embodiment comprising a dual-wavelength light source is provided for measuring 
relatively high levels of birefringence. 

Also provided is a technique for selecting the most accurate and efficient one 
of a number of approaches to determining birefringence properties of DUV optical 
elements, the selection depending upon the estimated level of the birefringence to be 
detected for a given optical sample. 

The birefringence properties (retardance magnitude and/or angular orientation) 
are precisely calculated. The system permits multiple measurements to be taken 
across the area of a sample to detect and graphically display variations in the 
retardance across the sample area. 

Other advantages and features of the present invention will become clear upon 
study of the following portion of this specification and drawings. 

Brief Description of Drawings 

Fig. 1 is a diagram of one preferred embodiment of the present invention 
showing a preferred arrangement of the optical components of a birefringence 
measurement system. 

Fig. 2 is a block diagram of the processing components of the system depicted 
in Fig. 1 . 



Fig. 3 is a diagram of another preferred embodiment of the present invention 
showing a preferred arrangement of the optical components of that birefringence 
measurement system. 

Fig. 4 is a block diagram of the processing components of the system depicted 
5 in Fig. 3. 

Fig. 5 is a graph depicting retardation curves for a sample measured at two 
different wavelengths in accord with on aspect of the present invention. 

Fig. 6 is a graph depicting retardation curves for a sample measured at three 
different wavelengths in accord with another aspect of the present invention. 
10 Fig. 7 is a drawing depicting a graphical display provided by the system of the 

present invention. 

Best Modes for Carrying Out the Invention 

One system for measuring birefringence is described in published 
International Patent Application No. W099/42796, as well as in US Patent 

15 application number 09/308,747, hereby incorporated by reference. That system 

employs a photoelastic modulator (PEM) for modulating polarized light that is then 
directed through a sample. The beam propagating from the sample is separated into 
two parts. These separate beam parts are then analyzed at different polarization 
directions, detected, and processed as distinct channels. The detection mechanisms 

20 associated with each channel detect the light intensity corresponding to each of the 

two parts of the beam. This information is employed in an algorithm for calculating a 
precise, unambiguous measure of the retardance induced by the sample as well as the 
angular orientation of birefringence relative to the fast axis of the sample. 
Considerations such as the nature of the light source required for retardance 

25 measurement at deep ultraviolet wavelengths (DUV) introduce the need for a 

somewhat different approach to birefringence measurement in the DUV environment. 

One preferred embodiment of the present invention uses a dual PEM setup to 
measure low-level linear birefringence in optical elements. This embodiment 
determines birefringence properties (both magnitude and angular orientation) that are 

30 the most important ones for CaF2 and fused silica suppliers to the semiconductor 



industry. This embodiment has specifically designed signal processing, a data 
collection scheme, and an algorithm for measuring low-level linear birefringence at 
very high sensitivity. 

As shown in Figure 1 , the dual-PEM setup 20 of this embodiment contains 
5 three modules. The top module comprises a light source 22, a polarizer 24 oriented 
at 45 degrees, and a PEM 26 oriented at 0 degrees. 

The bottom module includes a second PEM 28 that is set to a modulation 
frequency that is different from the modulation frequency of the first PEM 20. The 
second PEM 28 is oriented at 45 degrees. The bottom module also includes an 

10 analyzer 30 at 0 degrees and a detector 32. 

The middle module is a sample holder 34 that can be mounted on a computer- 
controlled X-Y stage to allow the scan of an optical element or sample 36. 

This embodiment (Figs. 1 and 2) employs as a light source 22 a polarized He- 
Ne laser at 632.8 nm. And, while the wavelength of this source is not DUV, the 

15 following is useful for explaining the general operation and analysis underlying the 
other dual-PEM embodiments explained below in connection with the DUV light 
sources that they employ. 

With continued reference to Fig. 1, the polarizer 24 and analyzer 30 are each a 
Glan-Thompson-type polarizer. A Si-photodiode detector 32 is used in this 

20 embodiment. Both PEMs 26, 28 are bar-shaped, fused silica models having two 

transducers. The transducers are attached to the fused silica optical element with soft 
bonding material. To minimize birefringence induced in the optical element, only the 
transducers are mounted to the PEM housing. The two PEMs 26, 28 have nominal 
resonant frequencies of 50 and 55 KHz, respectively. 

25 With reference to Fig. 2, the electronic signals generated at the detector 32 

contain both "AC" and "DC" signals and are processed differently. The AC signals 
are applied to two lock-in amplifiers 40, 42. Each lock-in amplifier, referenced at a 
PEM's fundamental modulation frequency (IF), demodulates the IF signal provided 
by the detector 32. In a preferred embodiment, the lock-in amplifier is an EG&G 

30 Model 7265. 
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The DC signal is recorded after the detector 32 signal passes through an 
analog-to-digital converter 44 and a low-pass electronic filter 46. The DC signal 
represents the average light intensity reaching the detector 32. As discussed next, the 
DC and AC signals need to be recorded at different PEM retardation settings. 

The theoretical analysis underlying the measurement of the birefringence 
properties of the sample 36 in this embodiment is based on a Mueller matrix analysis, 
and is discussed next for this dual PEM-single detector embodiment of Figs 1 and 2. 

For clarity, the Mueller matrices for three of the optical components in Fig. 1 
are shown below. The sample 36 in the optical arrangement, with a magnitude of 8 
and an angle of the fast axis at p, has the following form: 
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The Mueller matrices of the two PEMs, with the retardation axes oriented at p 
= 0° and 45° are, respectively: 
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where 81 and 82 are the time varying phase retardation of the first PEM 26 
and second PEM 28 (81 = 81 0 sina>it and 82 = 82 0 sinco 2 t; where coi and ©2 are the 
PEMs' modulating frequencies; 81 Q and 82 0 are the retardation amplitudes of the two 
20 PEMs). 

Using the Mueller matrices of the optical components in the set-up shown in 
Fig. 1, the light intensity reaching the detector 32 is obtained as follows: 




eqn. (1) 



where Io is the light intensity after the polarizer 24 and K is a constant that 
represents the transmission efficiency of the optical system after the polarizer. 

The functions of sin51 and cos51 in equation 1 can be expanded with the 
Bessel functions of the first kind: 



where k is either "0" or a positive integer, and J2k+i is the (2k+l) order of the 
Bessel function; and 



where J 0 is the 0* order of the Bessel function, and J2k is the (2k) order of the 
Bessel function. 

Similar expansions can be made for sin82 and cos52. 

Substituting the expansions of sin51, cos81, sin82 and cos52 into equation (1) 
and taking only up to the second order of the Bessel functions, we obtain the 
following terms: 



sin S\ = sin( Sl 0 sin(co x t)) = £ 2J 2k+l (Sl 0 ) sin((2A: + l)a> x t) 



eqn. (2) 



cos S\ = cos(£l 0 sin(^O) = J 0 (S\ 0 ) + £ 2J 2k (Sl 0 ) cos((2£)oV) 
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term (2) 
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term (4) 



The first parts of terms (3) and (4) can be used for determining linear 



5 retardance at low levels (below 7i/2 or a quarter- wave). Term (2) is useful for 
determining linear retardance at higher levels (up to n or a half-wave). Term (1) 
contains DC terms that relate to the average light intensity. 

The IF AC signals on the detector 32 can be determined using the lock-in 
amplifiers 40, 42 referenced at the PEMs' first harmonic (IF) frequencies. The lock- 
10 in amplifier will effectively exclude the contributions from all other harmonics. The 
IF signals measured by the lock-in amplifiers 40, 42 for the two PEMs 26, 28 are: 



where V2 results from the fact that the output of a lock-in amplifier measures 
the root-mean-square, not the signal amplitude. It is seen from eqn (4) that the 
15 maximum values of J 0 (51 0 )2Ji((52o) and Jo(52 0 )2Ji((51 0 ) will lead to optimal results 
for the output of the lock-in amplifiers. When the AC signals are collected, the 
retardation amplitudes of both PEMs are set to be 1.43 radians to optimize the AC 
signals. 



V2-F 11F = —^-J 0 (8\ 0 )- 2J X (82 0 ) cos(2 />) sin S 
KT 

V2-K 21/r =— <Ly 0 (,52 0 ). 2J, (<?l 01 )sin(2p)sin6 



eqn. (4) 



The DC signal can be derived from term (1) to be: 




eqn. (5) 
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where any term that varies as a function of the PEMs' modulation frequencies 
is omitted because they have no net contribution to the DC signal. The low-pass 
electronic filter 46 is used to eliminate such oscillations. 

Within small angle approximation (sinX = X and sin 2 X = 0 when X is small), 
5 V DC is independent of the sample's retardation and thus represents the average light 
intensity reaching the detector. However, when a sample with retardation above 30 
nm is measured, the V DC as shown in equation (5) will generally be affected by the 
magnitude and angle of the retardance. Thus, the measured DC signal will not be a 
true representation of the average light intensity. In this case, the most 
10 straightforward method is to set both Jo(51o) and Jo(52o) equal to "0". The DC signal 
then becomes: 

V nc=^J- eqn.(6) 

In this embodiment, the PEMs' retardation amplitude was selected as 
51o = 52 0 = 2.405 radians (0.3828 waves) for recording the DC signal. At such PEM 

15 settings, Jo(51o) = Jo(52 0 ) = 0. Therefore, the DC signal, independent of p or 8, truly 
indicates the average light intensity reaching the detector. 

As seen, this method requires recording AC and DC signals at different PEM 
settings and thus has a slower measurement speed (about 2 seconds per data point). 
This method affords high accuracy measurement of linear retardance above 30 nm. 

20 When speed is critical, an alternative method can be used. If the DC signal is 
collected at 81o - 52 0 = 01.43 radians, where the AC signals are recorded, the 
measured retardance of a sample, using the ratio of AC to DC, will depend on the 
sample's angular orientation. However, the DC term is well defined in equation (5). 
It is, therefore, possible to reduce the angular dependence of retardance by iteration of 

25 calculation for both retardation magnitude and angle. 

In order to eliminate the effect of light intensity variations due to light source 
fluctuations and the absorption, reflection and scattering from the sample and other 
optical components, the ratio of the IF AC signal to the DC signal are used. The 
ratios of AC signals to the DC signal for both PEMs are represented in equation (7): 



10 



^- = J 0 ( Sl 0 ) • 2J X (S2 0 ) sin S cos(2 p) 



eqn . (7) 

2 V V " F =JoWo)- 2^, (Sl 0 ) sin 5 sin(2 p ) 

V DC 

Defining Ri and R 2 as corrected ratios for both PEMs yields: 
: R, = sin 8 cos(2p) 

J 0 (Sl 0 )-2J x (S2 0 )-V DC " eqn (g) 

^ = R 2 =sin Jsin(2p) 

J 0 {52 0 ) 2J X {8\ 0 )V DC 

Finally, the magnitude and angular orientation of the birefringence are 
5 expressed as: 
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eqn. (9) 



where 8, represented in radians, is a scalar. When measured at a specific 
wavelength (i.e., 632.8 nm), it can be converted to retardation in nanometers: dnm = 
drad(632.8/(27i)). 

10 It should be emphasized that equations (9) are specifically developed for small 

linear birefringence due to the use of arcsine function in determining linear 
birefringence. Therefore, this method described here has a theoretical upper limit of 
7t/2 or 158.2 nm when using 632.8 nm laser as the light source. 

The signals at both PEMs' modulation frequencies depend on the orientation 

15 of the fast axis of the sample (see equation (6)), and the final retardation magnitudes 
are independent of the fast axis angles (see equation (9)). To achieve this angular 
independence of retardation magnitude, it is important to accurately orient all optical 
components in the system (as well as those of the embodiments described below). 

In this embodiment, the first PEM's optical axis is used as the reference angle 

20 ("0°"). All other optical components in the system are accurately aligned directly or 
indirectly with this reference angle. With the first PEM 26 being fixed, the following 
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procedures ensure the accurate alignment of all other optical components in the 
system: 

1 . With the second PEM 28 (50KHz) being turned off and the first PEM 26 

(55KHz) operating at quarter- wave peak retardation, the polarizer 24 
and analyzer 30 are approximately oriented at +45 degrees and -45 
degrees, respectively. 

2. Rotate the polarizer 24 in fine increments while monitoring the 2F (1 10 

kHz) signal from lock-in amplifier 40. When the 2F signal reaches its 
minimum (usually <0.05 mV with a lock-in sensitivity of 1 mV), read 
precisely the angle on the rotation stage of the polarizer 24. 

3. Rotate the polarizer 24 by precisely 45°, which is the correct position for 

the polarizer. 

4. Once the orientation of the polarizer 24 is correctly established, rotate the 

analyzer 30 in front of the detector 32 until the 2F (1 10 kHz) signal 
from lock-in amplifier 40 reaches its minimum. 

5. With the first PEM 26 (55KHz) being turned off and the second PEM 28 

(50KHz) operating at quarter- wave peak retardation, rotate the second 
PEM until the second 42 lock-in amplifier's 2F (100 kHz) signal 
reaches its minimum. 
When the optical components are misaligned, retardation magnitude shows 
specific patterns of angular dependence. 

The birefringence measurement of the present embodiment is specifically 
designed for accurately measuring low-level linear birefringence. In order to 
accurately measure such low levels of retardation, it is critical to correct for the 
existing residual linear birefringence of the instrument itself (instrument offset) even 
when high quality optical components are used. 

The instrument offset is primarily due to the small residual linear 
birefringence in the PEMs (on the order of 0.1 nm). To correct the system offset, an 
average of several measurements without any sample is first obtained. The 
instrument offsets are corrected in the software when a sample is measured. Notice 
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that such corrections should only be done when the ratios are calculated using 
equations (8), not on the final results of 8 and p, eqn. (9). The instrument offsets 
should be constants (within the instrumental noise level) unless there is a change in 
either the alignment of optical components or laboratory conditions such as 
5 temperature. It is prudent to check the instrument offsets with some regularity. 

This offset correction works within the limit of small retardance when the 
Mueller matrices of retardance commute. In practice, this is the only case where an 
offset correction is needed. Since the residual retardation in the PEMs is so small (on 
the order of 0.1 nm), offset correction will not be necessary when measuring 
10 retardation higher than 50 nm. 

The foregoing embodiment was specifically designed for measuring low-level 
retardance (up to a quarter- wave of the light source's wavelength, i.e. 158 nm for a 
633 nm He-Ne laser; 39 nm for the 157 nm light). 

The next described embodiment, illustrated in Figs. 3 and 4, is suitable for 
15 accurate measurements of relatively higher levels of retardance. This is important 
because a commonly used optical element in a DUV environment is CaF 2 , which has 
an intrinsic birefringence of about 12 nm/cm along one crystal axis. Thus, any such 
sample with a few em's thickness will produce retardation higher than the just- 
mentioned 39 nm limit, thereby requiring a system that can measure such relatively 
20 high values of retardation, which system is described next. 

As shown above, the prior-described embodiment uses the IF signal from one 
PEM 26 (Ml IF) and the IF signal from the other PEM 28 (M21F) to determine up to 
quarter-wave (e.g., 39 nm) retardation. By processing the electronic signal of the 
detector at the frequency of (Ml IF + M21F), the retardation range can be extended to 
25 half-wave of the light source's wavelength. 

The theory of this extension has already been provided in the earlier 
embodiment, which shows the development of terms (1) - (4) above. 

The important equations relating linear retardation to the detector signal at the 
30 frequency of (Ml 1F+M21F) are listed below: 
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5 The retardation noted in equations (10) is from 0 to half wave. The fast axis 

angle is determined with the IF data. Selectively combining the IF data and the 
Ml 1F+M21F data optimizes the determination of retardation from 0 to half- wave. 
For example, the Ml 1F+M21F data is used with the IF data for determining 
retardation around quarter- wave where the IF data is not accurate. Also, the IF data 

10 is used with the Ml 1F+M21F data to calculate the retardation around 0 or half-wave 
retardation where the Ml 1F+M21F data is not accurate. 

With reference to Fig. 3, the optical setup 120 for this embodiment is in many 
respects the same as that described in connection with the embodiment of Figure 1, 
including a polarizer 124 oriented at 45° and a PEM 126 at 0°. The system also 

15 includes a second PEM 128 that is set to a different modulation frequency (than the 
first PEM) and is oriented at 45 degrees, an analyzer 130 that is oriented at 0° and a 
detector 132. A sample holder 134 is mounted on a computer-controlled X-Y stage to 
allow the scan of a sample 36. Some differences in the structure and operation of 
these components, as compared with those of the earlier described embodiment, are 

20 described more fully below. 

Figure 4 shows the electronic signal processing block diagram of the present 
embodiment. 

Unlike the prior embodiment, the embodiment of Fig. 3 incorporates a light 
source 122 that is capable of generating beams of different wavelengths in the DUV 
25 region. These beams are collimated 123, and separately directed through the sample 
136 and processed as described more below. 
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Before describing the preferred and alternative mechanisms for providing the 
light source 122, the rationale underlying the need for a dual or multiple wavelength 
light source is discussed. 

With reference to Fig. 5, a system configured, as the system illustrated in Fig. 
5 1 , to operate at a single wavelength only gives correct and unambiguous retardation 
measurements at low levels; namely less than one-quarter wavelength. (Occasionally 
the wavelength symbol lambda (k) is hereafter used in lieu of the term wavelength.) 
It will be appreciated, however that without knowing in advance that the retardance 
value (magnitude) will be within the zero to quarter-wave range, an ambiguity will be 
10 present when the actual retardance value is calculated. 

In particular, the graph of Fig. 5 shows on its ordinate the measured retardance 
values (determined from the analysis presented earlier). The abscissa shows actual 
retardance levels. Considering a single wavelength system, the intensity-related 
signals provided to the computer 48 (Fig. 2) and based upon a single-wavelength light 
15 source will correspond to the single wavelength trace 50 on the graph of Fig. 5 

(ignoring for the moment the dashed-line second wavelength trace 52 described later). 

Considering this "M" shaped, single wavelength line 50 of Fig. 5, one can see 
that for a given measured (by the system) retardance value 8m, there is associated with 
that wavelength four possible actual values of retardance 5. This ambiguity appears 
20 on the abscissa at 5, (k/2 - 5), (7J2 + 5), and (A, - 5), within the range 0 to X (633 nm 
for this example). For example, when the instrument gives a reading of 40 nm for a 
measured sample, the actual retardation of the sample could be 40 nm, 277 nm, 357 
nm, or 593 nm. The actual retardation of the sample could also be the above values 
plus any multiple order of full waves, which cannot be resolved with a single- 
25 wavelength instrument. 

In accord with the embodiment of Fig. 3, two wavelengths (X\ and X 2 ) are used 
for the retardation measurement, and the possible measurements that are based on 
these wavelengths appear in the graph of Fig. 5 as solid line 50 and dashed line 52. 
In a preferred embodiment (Figs. 3 and 4) the light source 122 comprises a 
30 deuterium lamp combined with a monochromator. The lamp irradiates a wide range 
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of wavelengths. The monochromator selects the wavelength that is desired for the 
particular birefringence measurement application (such as 157 nm +/- 10 nm). It is 
contemplated that other lamps such as mercury lamps and xenon lamps can be used 
for birefringence measurements in different spectral regions. 

As shown in Fig. 5, the same sample has two distinct retardation curves 50, 52 
when measured at the two different wavelengths (solid line for X\; dashed line for X 2 ) 
provided by the light source 122 of this (Fig. 3) embodiment. The four points that 
reflect the ambiguity at 8, (X\/2 - 6), (k\/2 + 8), and (X\ - 8), which are labeled as 1,2, 
3, and 4 respectively, when measured with only wavelength X\, have separate 
retardation values (labeled as 1 2/ 3/ and 4') when measured at X 2 . 

To determine the correct actual retardance of the sample in this situation, the 
computer 148 of the present invention is programmed to carry out the following 
algorithm: 

1. Measure the sample at X\ 9 which gives (8 M , Pm); 

2. Measure the same sample at X 2 , which gives (S' M , P'm); 

3. Using the measurement results from the two wavelengths, the 
sample's actual retardation can be determined: 

3.1. If 8 M = S'm, the sample's actual retardation is 8 = 8 M = 8' M 
and p = pM= P'm; 

3.2. If S' M = 8 M + A/2 (A = X2 - XI and X2 > XI), the sample's 
actual retardation is 8 = Xl/2 - 8 M and p = p M ; 

3.3. If 8'm = Sm - A/2, the sample's actual retardation is 8 = Xl/2 + 
8 M and p = p M + 90 degrees; 

3.4. If 8'm = 8m + A, the sample's actual retardation is 8 = AT - 8m 
and p = p M + 90 degrees. 

As shown in Fig. 5, around the peaks, valleys and intersections of the two 
curves 50, 52 the conditions are different. When the wavelengths X\ and X 2 are 
selected to be sufficiently different, for example, X 2 being about 20% of the other, 
lower wavelength X\, measurement results at both wavelengths can be used to 
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determine unambiguously what the actual retardation is within 1 full wavelength of 
the retardation at the longer wavelength. 

One such case is when the actual retardation is close to X\/4 or 3X\/4, For this 
case, retardation measured with X\ is not accurate due to the use of the arcsine 
5 function in calculating the retardation. However, the measurement result with X 2 is 
accurate since the wavelengths are set far enough apart. The X 2 results can be used to 
determine the actual retardation. The significant difference of the X 2 results in the 
vicinity of X\/4 and 3X\/4 can be used to distinguish which is the actual retardation. A 
preferred procedure for use in this case is as follows: 
10 1 . Measure the sample at XI, which gives (8 M , Pm) where 5m is close to 

Xl/4 (say within 5%); 

2. Measure the same sample at X2, which gives (8'm, p'm); 

3. if 5'm is close to Xl/4 (say within 5%), the sample's actual 
retardation is close to Xl/4 — (8 = 8' M , p = p'm); 

15 4. if 8' M is close to Xl/4 - A/2 (say within 5%), the sample's actual 

retardation is close to 3X1/4. In this case, the sample's actual 
retardation is -- (8 = 8' M + A/2 + Xl/2; p = p' M + 90). 
Similarly, when the actual retardation is close to X 2 /4, 3X 2 /4, X\/2, X 2 /2, X\ 9 X 2 , 
or the intersections of the two curves 50, 52, the sample's actual retardation can be 
20 determined using the combination of measurement results from X\ and X 2 . 

In one embodiment, such as where the sample is under consideration is 
comprised of CaF 2 , X\ is selected to be 157 nm and X 2 may be, for example, 165 nm. 

It is noteworthy here that other light sources are contemplated in this 
embodiment. For example, the light source can be two or more separate lasers at 
25 different wavelengths. A switching device, such as a flip mirror, can allow the 

individual beams to pass to the sample, one at a time. Alternatively, the source can be 
a tunable laser that offers multiple wavelengths. Wavelength selection can then be 
determined by the computer-controlled system. 

Another choice of light source is a laser that emits multiple wavelengths 
30 simultaneously. To select the wavelengths for birefringence measurement, one can 
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employ an optical filter wheel to selectively pass the proper wavelengths. A filter 
wheel contains multiple optical filters mounted to the wheel. Rotation of the wheel 
allows a certain optical filter to be inserted into the path of the light beam. 

It is also possible to use a broadband light source combined with a filter wheel 
5 or wheels to select the desired wavelengths. Different types of optical filters, 

including high-pass, low-pass, and band pass filters, can be used in the filter wheel. A 
combination of filter wheels can be applied when necessary. 

The computer 148 is used to control and coordinate selecting wavelengths 
from the light source, as well as driving the PEMs at an optimal level for measuring 
10 birefringence and collecting data at an optimized sequence and calculating the final 
results. 

It is pointed out that although Fig. 4 shows two lock-in amplifiers 140, 142, 
that number can be different. For example, the use of one lock-in amplifier to detect 
sequentially the signals at different frequencies is also contemplated. Three lock-in 
15 amplifiers to detect Ml IF, M21F, and Ml 1F+M21F signals simultaneously are also 
contemplated. Once can also use a combination of sequential and simultaneous 
measurements. 

When lock-in amplifiers are used in the system (Fig. 4), the method described 
above only requires the use of IF data from both PEMs. When the (Ml IF + M21F) 

20 data is collected in addition to the IF data of the two PEMs, the range of measurable 
birefringence is extended to half wave at each wavelength used. In this case, the 
retardation curve at each wavelength becomes one triangle with its maximum at half- 
wave of the corresponding wavelength, as compared to two triangles for each 
wavelength with its maximum at quarter- wave of the corresponding wavelength as 

25 shown in Fig. 5. Collecting the extra data simplifies the analysis. 

When the birefringence dispersion at the chosen wavelengths is not negligible, 
this effect must be taken into account. For a specific sample, the birefringence at two 
chosen wavelengths is related by a constant determined by the material's dispersion. 
When the dispersion is unknown, it can be measured with a calibrated birefringence 

30 measurement system. Correction of birefringence dispersion is necessary for 
measurements around 157 nm. 
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In instances where the sample 136 has retardation higher than 1 full wave 
(A,l), the relationship described above is modified for determining the actual 
retardation. In this case, it is preferable that XI and X2 are close so that when the 
sample is measured at the two wavelengths, the retardation measured will be at the 
5 same order, i.e. mA,l + 51 and mk2 + 82. It is easy to select wavelengths when a 
white light source combined with a monochromator is used. 

It is possible that there is still ambiguity for some discrete data points when 
very large retardation values are involved, such as retardation satisfying (m+l)A,i = 
mX.2. In such cases, a third wavelength can be used to tell them apart. 

10 More particularly, when only two wavelengths are used, the two curves 

representing the two wavelengths intersect at two positions (Fig. 5). There are 
ambiguities in determining the actual retardation at those two positions. In Fig. 6, 
these two intersections are labeled I and II. The actual retardation for those two 
positions are 5b and 5d, respectively. The measured retardation for the two positions 

15 are 51 and 511. There are ambiguities between 5a and 5b as well as between 5c and 5d. 
One solution to resolve the ambiguities is to use a third wavelength, as shown by the 
retardance curve 54 in Fig. 6, which curve 54 is associates with the third wavelength. 
Using the third wavelength gives very different results, 81' and 5IF, for the positions I 
and II, respectively. 

20 It is important to point out that there are only a very limited number of data 

points where there is any such ambiguity. For example, two wavelengths are 
sufficient to measure all retardation below XI2 without any ambiguity. A third 
wavelength is only required when the measured retardation are at 51 and 511 (within 
the error of the measurement in reality) to determine the actual retardation. 

25 Upon further consideration of Fig. 6, it will be appreciated that, as the actual 

retardation increases, it is possible that all three curves intersect at the same position. 
In this case, XI, X2, and X3 will all give the same measurement result. A fourth 
wavelength could then be used to resolve the ambiguities and determine the actual 
retardation. Similarly, in rare cases, more than four wavelengths may be required for 

30 resolving ambiguities for very high levels of retardation. 
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The foregoing analyses provide various processes for measuring birefringence. 
Some processes or modes are more efficient than others, depending generally upon 
the magnitude of the retardation to be measured. What follows is a useful method for 
selecting the appropriate process for given ranges of estimated birefringence. 

1. Preferred process for measuring small linear birefringence (<0.35 radian or 

35 nm using 632.8 He-Ne laser) 

When the samples studied have small linear birefringence, the DC signal can 
be derived from term (1) to be: 

V DC = ^ jl + J 0 (Sl 0 ) • J 0 (S2 0 ) • sin(4/>) sin 2 (|J J eqn. (5) 

where any term that varies as a function of the PEMs' modulation frequencies 
is omitted because they have no net contribution to the DC signal. A low-pass 
electronic filter is used to eliminate such oscillations. Within small angle 
approximation (sinX = X and sin 2 X = 0 when X is small), V D c is independent of the 
sample's retardation and the PEMs' peak retardation setting, and thus represents the 
average light intensity reaching the detector. 

KT 

V D c=~f eqn * (6) - 
Therefore, the PEMs' peak retardation setting can be set at 1 .43 radians to 
maximize the values of Jo(51o)2Ji((52o) and Jo(52o)2Ji((51o), thus to optimize the IF 
AC signal recovery using lock-in amplifiers. At the PEMs' peak retardation setting of 
1.43 radians, the effect of the magnitude of a retardance in the worst case (sin(4p) =1) 
on V DC is such that the maximum error for the DC signal will not exceed 1% and 2% 
for retardation below 35 nm and 50 nm, respectively. 

When the linear retardance in a sample is small, an instrument based on 
detecting both the DC signal and IF AC signals at PEMs' peak retardation setting of 
1.43 radians, offers high sensitivity and fast measurements. The sensitivity is better 
than 0.005 nm using a 632.8 nm He-Ne laser as the light source. The speed of data 
collection is effectively determined by the time constant setting on the lock-in 
amplifiers. 
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2. Preferred process for measuring relatively higher linear retardation (0.35 
radians to n/2) 

When a 632.8 nm He-Ne laser is used as the light source, the effect on the DC 
signal from a sample with retardation above 35 nm should not be neglected. In this 
5 case, the most straightforward method is to set both Jo(51o) and Jo(52 0 ) equal to "0". 
At the PEMs' retardation amplitude 51 0 = S2 0 = 2.405 radians (0.3828 waves), the DC 
signal, independent of p or 5, truly indicates the average light intensity reaching the 
detector. 

However, it is clear from equation (4): 
KI 

V2 • V hlF = — °- J 0 (£1 0 ) • 2J X (S2 0 ) cos(2/?) sin 8 
10 2 eqn. (4) 

V2 - V 2 XF =^J 0 (S2 0 ) • 2J X (Sl 0l ) sin(2p) sin 8 



that the IF AC signals from both PEMs cannot be collected at PEM settings of Jo(51o) 
= Jo(52 0 ) = 0 at which the DC signal is recorded. For measuring higher levels of 
linear retardance (35 nm or 0.35 rad. to 158 nm or n/2), the method described here 
15 requires recording AC and DC signals at different PEM settings, thus has a slower 
measurement speed (about 2 seconds per data point). 

When speed is critical, an alternative method can be used. If the DC signal is 
collected at 81 0 = 52 0 = 1.43 radians where the AC signals are recorded, the measured 
retardance of a sample, using the ratio of AC to DC, will depend on the sample's 
20 angular orientation. However, the DC term is well defined in equation (4). It is, 
therefore, possible to reduce the angular dependence of retardance by iteration of 
calculation for both retardation magnitude and angle. 

3. Preferred process for measuring linear retardation up to n or 3 16.4 nm using 

632.8 He-Ne laser 

25 When a 632.8 nm He-Ne laser is used as the light source, a combination of the 

above mentioned processes and the use of term (2) can be employed to determine 
linear retardance up to half wave (3 16.4 nm) or 7t. This process will involve: 
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1. measuring retardation (both magnitude and angle of fast axis) up to 035 rad. 
using process 1 

2. measuring retardation (both magnitude and angle of fast axis) from 0.35 to n/A 
using process 2 

5 3. for retardation from n/4 to 37i/4, measuring retardation magnitude using term 

(2) and measuring retardation angle using process 2 
4. measuring retardation from 37i/4 to n using process 2 with modified algorithm 

(Sactual = 71 - 8 if) 

The combination of different processes will give the most accurate 
10 measurement results for linear birefringence measurement from 0 to n. 

In any of the preferred mode, the light source can be a variety of choices as 
proposed earlier. 

The measured retardance values can be handled in a number of ways. In a 
preferred embodiment the data collected from the multiple scans of a sample are 

15 stored in a data file and displayed as a plot on a computer display. One such plot 100 
is shown in Fig. 7. Each cell 102 in a grid of cells in the plot indicates a discrete 
location on the sample. The magnitude of the retardance is depicted by color coding. 
Here different shadings in the cells represent different colors. In Fig. 7, only a few 
different colors and cells are displayed for clarity. It will be appreciated, however, 

20 that a multitude of cells can be displayed. The legend 104 on the display correlates 
the colors (the color shading is omitted from the legend) to a selectable range of 
retardance values within which the particular measurement associated with a cell 102 
falls. A line 106 located in each cell 102 extends across the center of each cell and 
presents an unambiguous visual indication of the full physical range (-90° to +90°) of 

25 the orientation of the fast axis of the sample at each sampled location. Thus, the 
orientation of the fast axis and the retardance magnitude measurements are 
simultaneously, graphically displayed for each location. With such a complete, 
graphical display, an inexperienced operator user is less likely to make errors in 
analyzing the data that are presented. 
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In a preferred embodiment, the just described retardance measurements are 
displayed for each cell as soon as that cell's information is computed. As a result of 
this instantaneous display approach, the operator observes the retardance value of 
each cell, without the need to wait until the retardance values of all of the cells in the 
5 sample have been calculated. This is advantageous for maximizing throughput in 
instances where, for example, an operator is charged with rejecting a sample if the 
birefringence value of any part of the sample exceeds an established threshold. 

Also illustrated in Fig. 7 is a contour line placed there as an example of a 
contour line that follows a common measured range of retardation magnitude. For 
10 simplicity, only a single one of several contour lines is shown for the low-resolution 
plot of Fig. 7. 

It will be appreciated that any of a number of variations for displaying the 
measured data will suffice. It will also be apparent from Fig. 7 that the means for 
setting parameters of how the sample is scanned (scan boundaries, grid spacing 

15 sample thickness, etc.) and the resulting data are conveniently, interactively displayed. 
Another approach to graphically displaying the retardance magnitude and 
orientation information provided by the present system is to depict the retardance 
magnitude for a plurality of locations in a sample via corresponding areas on a three- 
dimensional contour map. The associated orientations are simultaneously shown as 

20 lines or colors in corresponding cells in a planar projection of the three dimensional 
map. 

While the present invention has been described in terms of preferred 
embodiments, it will be appreciated by one of ordinary skill in the art that 
modifications may be made without departing from the teachings and spirit of the 
25 foregoing. For example, it is also contemplated that other source wavelengths may be 
used with the embodiments described above. Representative ones of such 
wavelengths are 193 nm and 248 nm. 



